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Abstract. In this paper we present a three-dimensional numerical model for the radio emission of Magnetic Chemically 
Peculiar stars, on the hypothesis that energetic electrons emit by the gyrosynchrotron mechanism. For this class of radio stars, 
characterized by a mainly dipolar magnetic field whose axis is tilted with respect to the rotational axis, the geometry of the 
magnetosphere and its deformation due to the stellar rotation are determined. The radio emitting region is determined by 
the physical conditions of the magnetosphere and of the stellar wind. Free-free absorption by the thermal plasma trapped 
in the inner magnetosphere is also considered. Several free parameters are involved in the model, such as the size of the 
emitting region, the energy spectrum and the number density of the emitting electrons, and the characteristics of the plasma 
in the inner magnetosphere. By solving the equation of radiative transfer, along a path parallel to the line of sight, the radio 
brightness distribution and the total flux density as a function of stellar rotation are computed. As the model is applied to 
simulate the observed 5 GHz lightcurves of HD 37479 and HD 37017, several possible magnetosphere configurations are found. 
After simulations at other frequencies, in spite of the large number of parameters involved in the modeling, two solutions in 
the case of HD 37479 and only one solution in the case of HD 37017 match the observed spectral indices. The results of our 
simulations agree with the magnetically confined wind-shock model in a rotating magnetosphere. The X-ray emission from the 
inner magnetosphere is also computed, and found to be consistent with the observations. 

Key words. Stars: chemically peculiar - Stars: circumstellar matter - Stars: individual: HD 37479, HD 37017 - Stars: mag- 
netic field - Radio continuum: stars 



1 . Introduction some cases the mass loss rate (^M = 10 ^ — 10 A'Iqyt 

d and the outflow terminal speed (voc — 600 km s^^) was de- 

- - 'Magnetic Chemically Pecuuar (MCP) stars are commonly . . .cu . i inoT /-■ . o ir inm\ 
, . , , ,T^ ^/t ^r^A ^ ^ , rived (Shore et al. 1987, Groote & Hunger 1997). 

characterized by strong (B ~ lO'^ - 10 Gauss) and ^ ^ m^r. . i u j £ 

. .. . , , „ • r- , , r,,, , , About 25% or MCP stars also show evidence or non- 

periodically-variable surface magnetic nelds. The observed , , .• • ■ ■ ^ , irim t • i 

^ . . , , , ,,, thermal radio continuum emission (Drake et al. 1987, Linsky 

magnetic vanabihty is related to the global magnetic held ^ , ,„r.-i t . i ,r^nA\ t. u i. . j . .u 

^ , ■' , . , . . , , . , , et al. 1992, Leone et al. 1994). It has been suggested that the 

topology or these stars which is mainly dipolar, with the mag- 
netic axis tilted with respect to the rotational axis. The observed 
variability is a simple consequence of stellar rotation (Babcock 
1949). 



physical processes responsible for the observed radio emission 
are strictly related to the interaction between wind and mag- 
netic field (Linsky et al. 1992, Babel & Montmerle 1997a, here- 
after BM97). The flat spectral index and the observed degree of 
These stars show evidence of anisotropic winds, as indi- . , i • . , • , . j • . r 

, , . . circular polarization have been interpreted in terms of gyrosyn- 

chrotron emission from continuously injected mildly relativis- 
tic electrons trapped in the stellar magnetosphere. 



cated by spectral observations of UV lines (Shore et al. 1987, 
Shore & Brown 1990). Theoretical studies have demonstrated 



that a stellar wind in the presence of a dipolar magnetic field ~i j- • • r • 11 11 

, ^ , ^ , , • , , The radio emission of MCP stars is variable with the 

can freely flow only from the polar regions, where the mag 



same period as the magnetic field variability, as shown by 
Leone (1991) and Leone & Umana (1993) for HD 37479 and 
HD 37017, and by Lim et al. (1996) for HR 5624. In particular, 
for HD 37479 and HD 37017, the minimum of the radio emis- 
sion coincides approximately with the zero of the longitudinal 
magnetic field, whereas the maxima coincide with the extremes 
of the magnetic field curves (Leone & Umana 1993). The ob- 
Send offprint requests to: C. Trigilio served modulation suggests that the radio emission arises from 

Correspondence to: c.trigilio@ira.cnr.it a stable corotating magnetosphere. The temporal variability of 



netic field lines are open (Shore 1987). Therefore in MCP stars 
the wind forms two polar jets, whereas at the latitudes near 
the magnetic equator the wind is inhibited and the matter is 
trapped, forming the so called "dead zone". The presence of 
jets and circumstellar matter can explain the emission features 
observed in the UV spectra and Ha wings (Walborn 1974). In 
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Alfven Surface 




Fig. 1. Cross-section of the magnetosphere of magnetic chemically peculiar stars: the stellar dipolar field (short dashed curves) is modified by 
the stellar wind; inside the Alfven surface (long dashed curve), the magnetic pressure exceeds the kinetic pressure of the wind, and the field lines 
maintain a dipolar geometry. The largest closed field line defines the "inner" magnetosphere, which confines the stellar wind ("dead zone"). 
Here the two stellar wind streams from opposite hemispheres of the stellar surface colhde, leading to a shock that produces an enhancement 
of the temperature of the gas and eventually to X-ray emission. This gas cools and acciunulates in the magnetic equatorial plane, forming a 
torus-like cloud (the two filled circles close to the star). The open magnetic field lines just outside the inner magnetosphere produce current 
sheets (shaded areas), where the electrons are accelerated up to relativistic energies. They eventually can propagate back toward the stellar 
surface following the field lines of the "middle" magnetosphere and emit radio radiation by the gyrosynchrotron emission process (wave 
arrows). The hot matter in the inner magnetosphere can absorb the radio radiation. The field lines close to the magnetic poles of the star are 
open ("outer" magnetosphere), and tend to a radial topology outside the Alfven surface. Only the wind from the polar caps can escape from the 
magnetosphere. 



the radio flux is probably related to the change of the orienta- 
tion of the emitting region in the space, due to the misalignment 
of magnetic and rotational axes. 

In this paper we present a three-dimensional gyrosyn- 
chrotron model developed with the aim of investigating the na- 
ture of the rotational modulation of the radio emission. This 
kind of study can be used to test the physical scenario proposed 
to explain the origin of radio emission from MCP stars. 

2. The model 

To explain the radio continuum and the X-ray emission from 
young magnetic B stars and MCP stars, Andre et al. (1988) 
proposed a model characterized by the interaction between 
the dipolar magnetic field and the stellar wind. Close to the 



star, where the magnetic pressure dominates over the kinetic 
pressure of the wind, the magnetic topology is mainly dipo- 
lar. As the magnetic field strength decreases, the kinetic pres- 
sure equals the magnetic one, defining the Alfven surface, and 
the wind opens the magnetic field lines, generating a "cur- 
rent sheet" in the magnetic equatorial plane (Havnes & Goertz 
1984). Here the electrons of the wind can be accelerated up to 
relativistic energies (Usov & Melrose 1992); eventuaUy those 
electrons return to the star, along the field lines, thus defining 
the gyrosynchrotron emitting region. In the magnetosphere we 
identify three main zones: 

i) the "inner" magnetosphere or "dead zone", completely in- 
side the Alfven surface, defined by the largest closed mag- 
netic line; here the wind is confined; 



C. Trigilio et al.: 3D model for radio emission from magnetic CP stars 



3 



ii) the "middle" magnetosphere, where the radio emission oc- 
curs, holding all the open magnetic field lines that generate 
the current sheets; 

iii) the "outer" magnetosphere, where the wind flows freely 
from the magnetic polar caps, with a magnetic topology 
almost radial out of the Alfven surface. 



Linsky et al. (1992) located as possible radio emitting regions 
two tori, one in the northern, the other in the southern magnetic 
hemisphere, inside the "middle" magnetosphere, the higher fre- 
quency being generated closer to the star, where the magnetic 
field is stronger. 

The recent discovery of coherent emission at 20 cm from 
HD 124224 by Trigilio et al. (2000), explained in terms of 
Electron Cyclotron Maser Emission (ECME), seems to con- 
firm the hypothesis of relativistic electrons accelerated in the 
current sheets and returning toward the star, where they are re- 
flected back. In fact the ECME can occur after the magnetic 
mirroring when, in particular conditions, a loss cone anisotropy 
can develop (Melrose & Dulk 1982). 

The basic scenario of the magnetosphere of a MCP star is 
sketched in Fig. 1. If this scenario is correct, it should also re- 
produce the rotational modulation of the observed radio emis- 
sion. We outline the basic steps of the model in Appendix A.l, 
while the stellar reference frame and its rotation are presented 
in Appendix A.2 and A. 3 . The numerical computation of the ra- 
dio emission toward the Earth is presented in Appendix B. 1 . In 
the following, we analyze the numerical sampling of the emit- 
ting region, and the physical conditions of the magnetosphere. 



2. 1. Sampling 

The space surrounding the star is sampled in a 3-D cubic grid 
(see Appendix A.2) and all the physical quantities relevant for 
the numerical computation of the radio emission (such as mag- 
netic field, electron density,...) are evaluated in each grid point. 

Since the frequency of the gyrosynchrotron radiation is pro- 
portional to the gyration frequency of the electrons, which in 
turn is proportional to the magnetic field intensity, different fre- 
quencies are mainly emitted in different regions of the magne- 
tosphere at different distances from the stellar surface, because 
of the strong gradient of i? in a dipole. Therefore the modelled 
region is a cube whose overall size is a function of the fre- 
quency. This cube is then sampled into smaller cubes, whose 
size is chosen so that any physical parameter inside them can 
be assumed constant. 

The reliability of the model is obviously greater if the size 
of the cube element is smaller but, on the other hand, the com- 
putational time increases strongly. The compromise between 
sampling and computational time is chosen on the basis of the 
convergence of the model (the flux for instance) as a function 
of the sampling toward an asymptotic value. In particular, to 
model the 5 GHz radiation from the MCP stars, the overall lin- 
ear size of the grid is 25 — 30 steUar radii, divided into about 
80 — 100 cube elements. 
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Fig. 2. Thermal electron number density (top panel), magnetic field 
strength (medium panel) and equatorial radius of the Alfven sur- 
face (low panel) versus magnetic longitude They are derived with 
= 4: Rq, Prot = 1 day. Bp = 10* Gauss, /? = 60°, M = 
10"^ Mqyt-^ and v^o = 600 km s"\ see Sect. 2.2 for details. 

2.2. Location of tiie Aifven surface 

The equatorial radius of the Alfven surface (ja) may be esti- 
mated by equating the kinetic energy density and the mag- 
netic energy density ^ (Usov & Melrose 1992): 

Ew = ^ (1) 

The term includes both radial and rotational components of 
the kinetic energy of the wind: 

Ew = ^pv"^ -h ^pW^rf^ (2) 

where p and v are the density and the speed of the wind respec- 
tively, Lo is the angular velocity of the star and d the distance 
of the generic point of the magnetic equator from the rotational 
axis, given by: 



1 — sin j3 cos^ C 



where /? is the obliquity of the magnetic field, Q is the magnetic 
longitude (defined to be zero in the Une located by rotational 
and magnetic equator planes) and r the radial distance from 
the centre of the star. The gas density p in the "outer" region 
can be estimated by the continuity equation: 



M 



P = 



(3) 



where M is the mass loss rate. The wind speed w is a function 
of the distance from the star. For the sake of simplicity, we use 
the wind speed function derived by Castor & Simon (1983): 



v{r) 



1 



(4) 
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Fig. 3. Section of the magnetic shell which deUmits the emitting mag- 
netospheric volume. The parameter I determines the outer magnetic 

field line, and, consequently, the size of the radio emitting region. In 
the outer magnetosphere (wind zone) the density decreases with dis- 
tance. Close to the star, a cooler absorbing material is present, which 
we ascribe to a torus of constant density ne,th and temperature Tth. 
The plasma temperature and density in the inner magnetosphere (Tp 
and rip) are discussed in Sect. 2.4. The parameter I also defines also 
the width Apoi of the polar caps, from where the wind can flow freely. 



This equation is similar to the velocity law v{r) ~ 
(1 — theoretically predicted by Friend & Abbott 

(1986) for a radiatively driven stellar wind. However, we em- 
phasize that the above dependence of the stellar wind on the 
distance applies only for the "outer" and "middle" magneto- 
sphere, while we expect a different velocity and density distri- 
bution in the "dead zone" (see Sect. 2.4), as the close magnetic 
flux tubes strongly influence the dynamics of the trapped wind 
(BM97). 

To estimate the magnetic energy density in the equatorial 
plane of the "inner" magnetosphere, we use the expression of 
the strength for a dipolar field: 



B 



1 



(5) 



where Bp is the strength of the magnetic field at the pole of the 
star. 

Substituting Eqs.(2) and (5) in Eq.(l) we derive the equa- 
torial Alfven radius ta as a function of the magnetic longitude 
( searching the real roots. The equatorial Alfven radius (rA), 
the corresponding magnetic field strength (Ba) and the ther- 
mal electron number density (ne,A) are plotted in Fig. 2 as a 
function of (. They are derived for i?, = 4 B.q, rotational 
period Prot = 1 day, polar field Bp = 10'* Gauss, obliquity 
/3 = 60°, mass loss M = 10"^ Mq yr"! and v^o = 600 km 
s~* (Drake et al. 1987). Figure 2 shows how the stellar rotation 
significantly deforms the magnetosphere. 



2.3. Energetic electron distribution and location of 
the "middle" magnetosphere 

The present model assumes that electrons are accelerated in 
equatorial current sheets, drawn in Fig. 1 as shaded areas out- 
side the Alfven surface. The volume of the magnetosphere 
where the energetic electrons can propagate is defined in Fig. 3 
by thick solid lines. 

The model assumes that, in each point of the "middle" mag- 
netosphere, the emitting electrons are isotropically distributed 
in pitch angle and have a power law energy distribution 

7V(7) oc 7Ve(7 - 1)-' (6) 

where 7 is the Lorentz factor and is the total number den- 
sity of the energetic electrons. As low energy cutoff, we choose 
7 = 1.2, corresponding to 100 keV. This choice does not in- 
fluence the computation of the gyrosynchrotron emission, as 
discussed by Klein (1987). We choose the exponent S consid- 
ering that in the solar flares, where an acceleration process as 
in MCP stars is supposed, microwave and hard X-ray emission 
are assumed to be generated by the same population of non- 
thermal electrons. The spectrum of the hard X-ray photons is 
proportional to when the electron population is de- 

scribed by Eq.(6). Since the observed hard X-ray spectrum in 
solar flares is a power law with an exponent ranging from —3 
to -4 (Longair 1992), 6 should be in tiie range 2.5 - 3.5. At the 
maximum of the solar flares, i.e. close to the impulsive event, 
the spectrum is generally quite flat (5 « 2 -3), and eventu- 
ally it becomes steeper, as the more energetic electrons lose 
their energy faster. If in MCP stars the electron acceleration 
occurs in the current sheets outside the Alfven radius, then the 
non-thermal electrons injected toward the star do not lose much 
energy via microwave emission, as they quickly cross the emis- 
sion region (tens of seconds). Therefore they will maintain the 
original flat energetic spectrum. In the following, we will con- 
sider a spectrum distribution index S in the range 2-4. 

The electron number density remains almost constant as 
the electrons propagate inside the magnetosphere, as a conse- 
quence of two combined effects (Hasegawa & Sato 1989): 1) 
particle losses due to magnetic mirroring and 2) compression 
due to the decreasing volume of the magnetic shell. 

We assume that the number density is a fraction of the 
thermal plasma density ne,A at the Alfven surface. This ratio 
is a free parameter in our model, since it depends on the effi- 
ciency of the acceleration process, and we assume A?e/ne,A in 
the range 10~''-1. 

The "middle" magnetosphere is defined by two magnetic 
lines (Fig. 3): the first is the largest closed line of the "inner" 
magnetosphere, tangent to the Alfven surface in the plane of the 
magnetic equator; the latter is defined by the size of the current 
sheet, which depends on the particular physical conditions of 
the magnetosphere, and which we consider a free parameter of 
our model. 

In a magnetic dipole, the equation of a field line is given by 
r = L cos^ A 

where L is the distance between the center of the star and the 
point where the field line crosses the magnetic equatorial plane; 
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Fig. 4. Effective magnetic field variations (Bohlender 1987) and 5 
GHz lightcurves of HD 37479 and HD 37017: (x) data from Drake 
et al. (1987), (•) Leone & Umana (1993), (o) this paper. 



A is the angle between this plane and the radius vector r. The 
middle magnetosphere is therefore identified by the points be- 
longing to magnetic field Unes having an equatorial distance 
L between ta and rx + l, where I is a free parameter in our 
calculations. 

2.4. Thermal plasma In the Inner magnetosphere 

The presence of thermal plasma trapped in the dead zone may 
affect the stellar radio emission because of the free-free ab- 
sorption (Andre et al. 1988). If its density is high enough, this 
plasma can partially absorb the radiation and this could result in 
an additional modulation of the radio lightcurve, since this ef- 
fect is greater for a particular orientation of the magnetosphere. 
Moreover, if the density of the plasma is higher close to the star. 



the highest frequencies, which are emitted in the inner part of 
the magnetosphere, will be strongly absorbed. In contrast, the 
lowest frequencies, generated far from the star, are less affected 
by the plasma absorption. Therefore, the density distribution 
inside the irmer magnetosphere can affect either the modula- 
tion of the radio Ughtcurve, or the spectrum of the source. In 
particular, it is important to define the extent of the absorbing 
matter into the inner magnetosphere, the density distribution 
and the temperature. 

Effects of the trapped plasma have been studied at other 
wavelengths, as it causes also the absorption of resonance lines 
and of continuum when the star is seen with the equator edge- 
on. Smith & Groote (2001), from the UV spectra of several 
magnetic early-B stars, found some absorbing material, which 
they interpreted as a "cloud" along the line of sight, character- 
ized by a column density (ricoi) in the range 10^2-23 cm-s^ 
temperature in the range 15 000 -31 000 K and average den- 
sity 10^^ ~ cm~^. 

With the aim of explaining the X-ray emission from IQ Aur, 
BM97 developed a magnetically confined wind-shock model, 
where the radiatively-driven winds from the two magnetic 
hemispheres collide, producing a shock that leads to an en- 
hancement of the temperature up to 10^ ^ ^ K in the post-shock 
region, where X-ray radiation is emitted. They also predicted, 
from their X-ray model, the existence of a "disk" along the 
equator, made of post-shock cooling material, which could in 
particular explain the 15.4'' rotational modulation of 0^ Ori C 
seen by ROSAT(^cih&\ & Montmerle 1997b). In this model, the 
cloud observed by Smith & Groote (2001) could be the inner 
part of the equatorial disk. Also Donati et al. (2001) interpreted 
the X-ray emission and the variabiUty of /3 Cep in the frame- 
work of this model. The physical scenario of such a magne- 
tosphere has been summarized by Montmerle (2001) (see also 
Fig. 1). 

Once a steady state of this process has been reached, the 
dynamics of the wind and the physical characteristics of the 
matter along lines of force of the magnetic field are ruled by 
the equilibrium of the total pressure. This implies that 



(7) 



with rip and Tp respectively the number density and the temper- 
ature of the plasma in the post-shock region and Pram = pv^ 
is the ram pressure of the wind. The density p and the speed 
V come from Eqs.(3) and (4). If M = 10~^ MQyr~^ and 
= 600 km s~^, then f>ram ~ 10~^ - 10~^ gcm~^s~^. 
Assuming that the radio absorbing region is the post-shock re- 
gion of BM97, the previous relations give ripTp k, 10^^ (in 
c.g.s.); if Tp « 10*^ K, rip « 10^ cm-^. 

If, in addition, the post-shock region were also responsible 
for the absorption of the UV lines and continuum (Smith & 
Groote, 2001), a column density of 10^^ cm~^ should lead to 
a characteristic size of the post-shock region given by ricoi/np 
that, considering the Eq. (7), is « lO^^/rip ~ lO^Tp. Therefore, 
to be responsible also for the absorption of both UV features 
and radio radiation, the X-ray emitting region (with Tp w 10^ 
K) should be very large (> 10^^ cm), at least some hundreds 
times the stellar radius or some tens the Alfven radius. In few 
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words, so high a column density cannot be reached inside the 
Alfven radius unless the density is very high. 

The other possibility is that the X-ray emitting region 
(BM97) is not responsible for the UV absorptions, and the only 
valid relation between Tp and rip is Eq.(7). If Tp « 10*^ " K, 
we find rip w 10^~^cm~^. Briefly, starting from the stellar 
surface, we have: 

a) a dense region (rie, th ~ 10^^ cm~^) at stellar temperature, 
(Tth « Teff) with a size of about 10^^ cm, less than the 
stellar radius: this region is responsible for the UV absorp- 
tion; for simphcity, a torus shape has been assumed to ac- 
count for the two equatorial matter distributions, with den- 
sity rie, th and temperature Tth assumed constant; 

b) a less dense region (rip « 10 ~ ^ cm""^) at high tempera- 
ture (Tp » 10^ ~ K), whose maximum extent is defined 
by the magnetic confinement (B"^ /Sir ss n^kT^). Since 
the magnetic field is given by Eq.(5), the maximum size 
is given by 



I 167rnpfcBrp j 

which gives L = 23 - 18 i?* respectively for = 10 000 
- 5 000 Gauss. In any case L is greater than the Alfven ra- 
dius (Fig. 2, lower panel). We can therefore assume that this 
region fills all the inner magnetosphere. 

Figure 3 also outlines this scenario. The actual situation can, 
however, be a little different since the above picture is vahd 
in the hypothesis of a non-rotating star. In fact, BM97 showed 
that the temperature of the plasma in the irmer magnetosphere 
remains almost constant if the star does not rotate, but can in- 
crease Unearly outward in the case of a fast rotation (with the 
dipole axis aligned with the rotational one, i.e with fi = 0). 
Since the plasma is in pressure equihbrium with the wind ram 
pressure, we have to add to the first term of Eq.(7) a kinetic 
term that accounts for the rotational energy, in a similar way as 
in Eq.(2). In any case, if the rotation is considered, the density 
should decrease outward. In the case that (3^0 the situation is 
more complex, but when /3 « 90° the rotational axis lies close 
to the magnetic equatorial plane, and the regions close to the 
rotational axis can be considered to be not rotating, while other 
regions can be rotating fast; the result is intermediate between 
a non-rotating and a rotating magnetosphere. 

The presence of thermal plasma has been considered in the 
calculation of the emitted radiation (see Appendix B.l), includ- 
ing free-free thermal emission and absorption (Dulk 1985) as 
weU as the Razin effect. 

3. Application to HD 37479 and HD37017 

The model developed in this paper has been used to simulate 
the radio lightcurves of two MCP stars: HD 37479 (=a Ori E) 
and HD 37017. 

The 5 GHz emission of these two stars appears to vary with 
the rotational period. Fig. 4 shows the effective magnetic field 
Be and the observed 5 GHz flux densities S5 as reported by 
Leone & Umana (1993). We add to these data the new re- 
sults of two VLA observations: HD 37479 was observed on 29 



Table 1. Stellar parameters 



Measured parameters 


D 


-Prot 


Bp i° 




[pc] 


[day] 


[Gauss] 




HD 37479 352^ 


1.19081'' 


6800'' 72* 


56* 


HD 37017 373^ 


0.901195" 


7700" 25* 


65* 


Derived parameters 


rA (min-max) 


rie.A (min-max) 








[xlO*^ cm"^] 




HD37479 13.18-15.75 


2.6-3.7 




HD 37017 12.59-16.36 


2.4-4.1 





' Hipparcos main catalogue 

" Bohlender et al. (1987) 
* Shore & Brown (1990) 



April 1995 (21:29-21:47 UT) and HD 37017 on 29 May 1995 
(18:42-19:06 UT). We reduced and analysed flie data foUowing 
the procedures used by Leone & Umana (1993). We measured 
as flux densities respectively 3.0 ±0.1 mJy at phase $ = 0.545 
and 1.9 ± 0.2 mJy at i> = 0.16, with rotational phases com- 
puted by using the ephemeris determined by Bohlender et al. 
(1987). Those data are reported in Fig. 4 as diamonds sym- 
bols. Leone & Umana (1993) noted that, while the radio data 
of HD 37017 can be fitted by a single sinusoidal like curve, 
for HD 37479 the radio variation can have a double-wave be- 
haviour, as for the variations of the helium abundance in those 
stars. They argued that the radio emission is probably related 
to \Be\, with the minimum corresponding approximately with 
the zero of the longitudinal magnetic field. 

Regarding the X-ray emission, that could arise from the 
inner magnetosphere, both stars have been observed in the 
ROSAT all-sky survey (Berghofer et al. 1996). While HD 37017 
is reported as a non-detection (IgLx < 30.7), HD 37479 is a 
doubtful detection (Drake et al. 1994) since it is a visual bi- 
nary with a hotter primary component closer than 1', Ukely the 
predominant contributor to the observed X-ray emission. 

To simulate the 5 GHz radio emission of these two MCP 
stars, we fix a stellar radius of i?* = 4 Rq (Shore & Brown 

1990) and an effective temperature Tdi = 24 000 K (Leone 

1991) . The other physical stellar parameters are reported in 
Table 1. We set the mass loss rate and the terminal wind speed 
respectively to M = 10^^ AfQyr^^ and v^c = 600 km s~^, 
according to Drake et al. (1987). The above mass loss rates 
are not the actual values, but refer to the case of a spherical 
wind. In fact, the wind is confined by the magnetosphere, lead- 
ing to the formation of the torus as discussed in Sect. 2.4, and it 
can escape only from the two polar caps defined by the largest 
closed field line (see Sect. 2.2 and Fig. 3). Using these observ- 
ing constraints we can univocally define the equatorial radius 
of the Alfven surface and the local thermal plasma densities. 
The maximum and minimum equatorial Alfven radius (rA) 
and the corresponding thermal electron number density (ne,A) 
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Fig. 5. HD 37017: examples of search of solutions. Light curve characteristics: maximum flux density (top panels) and ratio A (lower panels) 
as a function of the magnetic shell thickness I normalized to the stellar Alfven radius ta ■ In these examples, the spectral index of the relativistic 
electrons {S) has been fixed to 2, and the temperature of the plasma at the base of the inner magnetosphere (Tp) to 10^ K. The rotation of the 
star has been considered. From left to right: solutions for increasing values of the plasma density (rip) at the base of the inner magnetosphere. In 
each pair of panels, Smax and A are computed for different number density of the emitting electrons (iVr) as percentage of the number density 
of thermal plasma at the Alfven point rio.A (different line styles). A valid solution is given if the computed S'max and A are within the observed 
values (shaded areas) for the same value of IjrK- The good solution is indicated by a vertical line and two squares. 



for HD 37479 and HD 37017 are also reported in Table 1. The 

actual mass loss rate Mact is computed considering that the 
polar caps have an aperture Apoi defined by the intersection 
between the stellar surface and the magnetic field line that 
identifies the "outer" magnetosphere (see Fig. 3). Using the 
above stellar wind parameters, we get Apoi « 15°, so that the 
area of the two caps is only 2% of the total surface, giving 
Mact ~ lO-"M0yr-i. 

The aim of our analysis is to derive possible configurations 
of the magnetosphere, represented by the free parameters of 
the model (Z, A^e> T-^ and rip), which result in simulated 
lightcurves at 5 GHz able to reproduce the observations. The 
search of these combinations of the free parameters, computing 
the radio flux at different values of the rotational phase, could 
be very time-consuming so, to restrict the computing time, we 
simulate the flux density of the two stars only for particular 
stellar orientations, coinciding with the extremes and the nulls 
of the magnetic field curve. The characteristic of the observed 
radio lightcurves may be summarized by the maximum flux 
density ^max and the ratio A = S'min/'S'max, that expresses the 
depth of the modulation, being minimum for high modulation. 

The lightcurve characteristics derived by fitting the observ- 
ing data for HD 37479 and HD 37017 are, respectively: 
Smax = 4.0 ± 0.2 mJy, A = 0.65 ± 0.06 and 
S'max = 2.6 ± 0.1 mJy, A = 0.58 ± 0.04. 

The simulations for the two stars are performed by using 
different combinations of the free parameters of the model: 



- I: equatorial thickness of the magnetic shell (Fig. 3), such 
thatZ/rA = 0.025-1; 

- A^e^ total number density of the non-thermal electrons, in 
the range 10^ - lO'^cm-^, wifli < ne,A; 

- 5: spectral index of the non-thermal electron energy 
distribution (Eq.(6)), in the range 2-4; 

- Tp and rip: temperature and number density of the plasma 
in the inner magnetosphere (the post-shock region of 
BM97); they are chosen so that TpUp = pram/A:B is in the 
range lO^"^ - 10^^ (in c.g.s.), with Tp = 10^ - lO'^ and 
rip = lO'' - 10^°cm-3; 

- rotation: if the star does not rotate, the thermal plasma 
density is considered constant within the inner magneto- 
sphere; if, on the contrary, we consider the rotation, the 
thermal plasma density decreases hnearly outward, while 
the temperature increases. Tp and rip are considered as the 
values at r = i?*. 

The behaviour of 5i„ax and A versus the magnetic shell 
thickness I is investigated for models computed assuming dif- 
ferent values of A^e. Tp, rip, and rotation. Fig. 5 shows some 
examples of the adopted method of searching for good com- 
binations of the free parameters. The shaded areas represent 
the observed values of S'max and A derived from the observed 
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Table 2. With and without rotational effects, the 5 GHz radio lightcurves of HD 37479 and HD 37017 are reproduced for the listed values of: 
spectral index of non-thermal electron energy distribution (5); temperature (Tp) and thermal electron density (np) of the plasma in the post- 
shock region as extrapolated to the stellar surface; equatorial thickness in Alfven radius units (/ /a a); and total number density of non-thermal 
electrons (A^ ). The asterisk in the last column indicates the possible solutions that reproduce the ahnost flat observed spectral index. 



Star 


Rotation 


(5 


Tp[K] 


rip [cm 


l/r-A 


jVr[cm 


Spectral index 


HD 37479 


No 


2 


1.0 


X 


lo- 


> 1.0 


X 


10'^ 


0.74 


6.3 


X 


10^ 


-0.40-0.10 






Yes 


2 


1.0 


X 


10^ 


1.0 


X 


lO'^ 


0.05 


3.1 


X 


10^ 


-0.18-0.31 






Yes 


2 


1.0 


X 


10^ 


2.0 


X 


10» 


0.40 


6.3 


X 


10^ 


-0.10-0.34 


* 




Yes 


2 


1.0 


X 


10^ 


3.0 


X 


10^ 


0.83 


4.5 


X 


10^ 


-0.07 - 0.21 


* 




Yes 


3 


1.0 


X 


10^ 


1.0 


X 


lO'' 


0.11 


3.1 


X 


10* 


-0.07-0.61 






Yes 


3 


1.0 


X 


10^ 


2.0 


X 


10« 


0.77 


9.4 


X 


10^ 


-0.04 - 0.56 






Yes 


4 


1.0 


X 


10^ 


1.0 


X 


10 


0.27 


3.1 


X 


10^ 


0.05 - 0.74 






Yes 


4 


1.0 


X 


10^ 


2.0 


X 


10» 


1.05 


1.5 


X 


10^ 


0.11-0.69 






No 


2 


1.0 


X 


10" 


5.0 


X 


10® 


0.36 


9.4 


X 


10^ 


0.12-0.75 






No 


2 


1.0 


X 


10" 


1.0 


X 


10^ 


0.74 


6.3 


X 


10^ 


0.09 - 0.54 






No 


3 


1.0 


X 


10® 


3.0 


X 


10® 


0.18 


3.1 


X 


10* 


0.06-0.87 






No 


3 


1.0 


X 


10® 


4.0 


X 


10® 


0.56 


1.5 


X 


10* 


0.09-0.92 






No 


3 


1.0 


X 


10® 


5.0 


X 


10® 


0.91 


1.2 


X 


10* 


0.12-0.99 






No 


4 


1.0 


X 


10" 


1.0 


X 


10® 


0.12 


7.8 


X 


10-' 


0.02-0.81 






No 


4 


1.0 


X 


10® 


2.0 


X 


10® 


0.12 


7.8 


X 


lO'^ 


0.11-0.93 






No 


4 


1.0 


X 


10® 


3.0 


X 


10® 


0.45 


3.1 


X 


10" 


0.16-0.96 






No 


4 


1.0 


X 


10® 


4.0 


X 


10® 


0.83 


2.2 


X 


10" 


0.21 - 1.05 






Yes 


4 


1.0 


X 


10" 


1.0 


X 


10® 


0.05 


1.5 


X 


10" 


-0.09 - 0.58 






Yes 


4 


1.0 


X 


10® 


1.0 


X 


10^ 


0.05 


1.5 


X 


10® 


-0.09 - 0.60 






No 


4 


1.0 


X 


10^ 


1.0 


X 


10^ 


0.07 


1.5 


X 


10® 


-0.04-0.71 






No 


4 


1.0 


X 


10^ 


1.0 


X 


10® 


0.07 


1.5 


X 


10" 


-0.04-0.71 






Yes 


4 


1.0 


X 


10^ 


> 1.0 


X 


10^ 


0.09 


7.8 


X 


10" 


-0.11-0.66 




HD 37017 


No 


2 


1.0 X lO" 


> 1.0 X 10^ 


0.59 


4.5 X 10^ 


-0.44 - 0.09 






Yes 


2 


1.0 X 10^ 


3.0 X lO** 


0.53 


3.1 X 10^ 


0.00-0.34 


* 




Yes 


3 


1.0 X 10^ 


2.0 X 10^ 


0.62 


4.7 X 10^ 


0.14-0.68 






Yes 


4 


1.0 X 10^ 


2.0 X 10^ 


0.67 


9.0 X 10* 


0.08-0.86 






No 


2 


1.0 X 10® 


5.0 X 10® 


0.05 


3.1 X 10^ 


0.27-1.10 






No 


2 


1.0 X 10® 


1.0 X lO'' 


0.59 


4.5 X 10^ 


0.22-0.60 






No 


3 


1.0 X 10® 


4.0 X 10* 


0.40 


9.4 X 10^ 


0.29-1.18 






No 


4 


1.0 X 10® 


3.0 X 10® 


0.36 


1.5 X 10^ 


0.29-1.12 





5 GHz lightcurves. For each value of N^, the good solutions 
must show the observed values of S'max and A for the same 
value ofl/vr. They are reported in Table 2 for the two stars. 



4. Discussion 

The model presented in this paper is used to test the scenario 
proposed to explain the origin of radio emission from MCP 
stars and, if possible, to derive quantitative information about 
the physical conditions of their magnetospheres. 



simulated spectra must be flat (say in the range -0. 1 - 0.4) as 
reported by Drake et al. (1987) and Leone et al. (1994). The 
solutions with the required spectral index are indicated by an 
asterisk in the last column of Table 2. 
The model is not computed at 1.4 and 22 GHz because: 

1) at low frequency the emitting region extends far from the 
star and close to the Alfven surface, where the geometry of the 
magnetosphere is not yet well known, and 

2) at high frequency the magnetosphere must be so closely 
sampled that computational times are prohibitive. 



First, the characteristic quantities 5i„ax and A of the 
lightcurves at 5 GHz of HD 37479 and HD 37017 are com- 
puted for different combinations of /, N^, 6, Tp, rip, including 
also the different behaviour of the thermal plasma in the irmer 
magnetosphere as the star rotates or not. AH the possible 
solutions are shown in Table 2. 

Second, for any possible solution, we compute the emerg- 
ing radio flux also at 8.4 and 15 GHz; the spectral index of the 



Third, we compute the radio lightcurves at 8.4 and 15 GHz 
with the aim of predicting the behaviour of our model, which 
should be tested in the future with multi-frequency observa- 
tions. 
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4.2. The effect of the electron energy spectrum 



10 



i 10 



HD 37479 



HD 37017 



1 10 
p [GHz] 

Fig. 8. Top panel: Observed spectrum of HD 37479 by Leone et al. 
(1994) (•) and its linear 2-15 GHz fit (dotted lines) are compared 
with the computed spectra at minimum and maximum radio emission 
for the two solutions marked with asterisks in Table 2. In the order they 
are listed, dashed and dashed-dotted lines represent these solutions. 
Bottom panel: with the previous symbols, the observed spectrum of 
HD 37017 is compared with the only solution fotmd. 



4. 1. Effect of the rotation 

The radio lightcurves at 5 GHz for the two stars can be repro- 
duced by several combinations of free parameters. Particularly, 
the effects of the rotation seem to be important. 

Neglecting the effects of rotation, i.e. assuming that the 
density and the temperature of the inner magnetosphere are 
constant, we find that Tp = 10^ K, the electron energy spec- 
trum must be hard (S = 2) and a small relativistic electron num- 
ber iVe is required to reproduce the observed radio emission. In 
this case rip must be quite high (> 1 x 10^ cm~^), in order to 
have the observed modulation of the radio lightcurve. Also if 
Tp = 10^ K we get solutions for both stars, for S = 2,3, 4; 
A^e increases as the electron energy spectrum becomes softer 
(S = 4), reaching almost the number density of the Alfven 
point (i.e. « 10® cm~^). rip is lower than in the previous case, 
being in the range 10® - 10^ cm~^. 

For Tp = lO'^ K, only the behaviour of HD 37479 can be re- 
produced, but only for J = 4 and very high values of N^. 

The situation is different if the rotation has been taken into 
account, considering that Tp increases linearly outward and rip 
decreases according to the relation Tprip = const. In Table 2 
their values extrapolated at i? = i?* are reported. 
For Tp = lO'' and S = 2,3, 4, the lightcurves of both stars are 
well reproduced. The number density at the base of the inner 
magnetosphere is about 10^ cm~^. 

For Tp = 10® K we get solutions only for HD 37479, only with 



5 = 4 but with 



ne,A- 



Simulations at Tp = lO'^ K do not show any possibility for 
HD 37017, while for HD 37479 they are limited to J = 4 and 
very high A''e. 



As already noted, a hard non-thermal electron energy popula- 
tion, such as 5 = 2, is able to reproduce the observed radio 
emission with small A^e (~ 10^ - 10^ cm~^). A softer popula- 
tion, say 5 = 4, even if able to reproduce the emission, requires 
a very efficient acceleration process, as almost all the electrons 
flowing with the stellar wind should be accelerated up to rela- 
tivistic energies. 

Our results indicate that a hard electron energy population is 
more probable. 

4.3. Simulation of radio lightcurves and spectra 

In order to get the spectral index, radio fluxes at 8.4 and 
15 GHz are computed for any solution at maximum and min- 
imum hghtcurve. The range of variation of the spectral index 
is also reported in Table 2. Only two solutions for HD 37479 
and one for HD 37017 give spectral indices similar to the ob- 
served ones. They have in all cases a quite hard electron energy 
distribution {S = 2) and are all consistent with a rotating mag- 
netosphere, where the density of the absorbing matter in the 
inner magnetosphere decreases outward quite hnearly, while 
the temperature increases. Radio spectra at maximum and min- 
imum flux are shown in Fig. 8 for the accepted solutions, to- 
gether with the spectral data from Leone et al. (1994). There is 
a good agreement, and this result encourages us to extend the 
model at lower and higher frequency. 

The simulated radio lightcurves for HD 37479 and 
HD 37017, corresponding to the three accepted solutions, are 
shown in Fig. 9. In the left panels, the 5 GHz simulations are 
superimposed on the measurements from the literature. For the 
first star, there is no significant difference between the two sim- 
ulated radio lightcurves and it is, therefore, impossible to dis- 
criminate between them by using only single frequency data. 
Even though radio Hghtcurves are not available for those two 
stars at other frequencies, we computed the synthetic curves at 
8.4 and 15 GHz in order to give a prediction of the behaviour 
in a wide range of frequency. The comparison between the re- 
sults of our model computed with multifrequency observations 
is therefore crucial. 

Both stars are unfortunately too far away (Table 1) to give 
us the opportunity to resolve the radio sources by using Very 
Long Baseline Interferometry (VLBI). In fact, at 5 GHz, the 
extent of the radio emitting region (Fig. 6) is about 10 stellar 
radii, corresponding to an angular size of about 1 .5 milliarcsec- 
onds (mas). The angular resolution of the VLBI arrays (with 
a maximum baseline of about 8000 km) is 1.5 mas, and the 
source cannot resolved. This is in agreement with the results of 
Phillips & Lestrade et al. (1988) who did not resolve HD 37479 
with high-sensitivity intercontinental VLBI. At 15 GHz the res- 
olution of the VLBI is 0.5 mas and the sources can be only 
barely resolved if observed at a rotational phase corresponding 
to the zero of longitudinal magnetic field B^, when the radio 
source has its maximum angular extent (a Uttle less than 10 
stellar radii). A closer MCP star should be observed in VLBI 
in order to test the morphology predicted by our model. 
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Fig. 9. Left panels: comparison of the observed 5 GHz data of Fig. 4 and the computed lightcurves corresponding with the accepted solutions 
of Table 2. Central and right panels: lightcurves at 8.4 and 15 GHz as predicted by the model. See Fig. 4 for the meaning of the symbols. 



4.4. X-ray emission 

The physical conditions of the inner magnetosphere play an im- 
portant role in the emerging radio flux because of the free-free 
absorption by plasma. The inner magnetosphere is modelled 
following the magnetically wind-shock model by BM97. By 
performing more than 56 000 simulations per star, changing all 
the free parameters of our model, we find only 3 solutions satis- 
fying all the observational constraints. All 3 solutions indicate 
that the effects of rotation are not negligible. For the two stars, 
we find that the temperature increases outward from about 10^ 
to 1.4 X 10^ K, while the density decreases from 2 - 3 x 10^ to 
2 X 10*^ cm^"*. As a by-product of our model, we can estimate 
the X-ray emission of the inner magnetosphere once its config- 
uration has been found according to the radio emission model. 
Following a procedure similar to that used for the computation 
of the radio flux, we compute for each element of the magneto- 
sphere the temperature and the number density in order to get 
the Bremsstrahlung emission coefficients for a thermal plasma. 
The emitted power is then integrated in the range of energy 
0.1 — 2.0 KeV, and the resulting X-ray luminosities are eval- 
uated as IgLx = 30.06 for HD 37479 and IgLx = 30.42 for 
HD 37017. These values are consistent with the hmits of detec- 



tion of the ROSAT all-sky survey, as in Sect. 3. However, with 
the new X-ray telescopes now available (Chandra and XMM), 
both stars should be easily detectable. In fact, the correspond- 
ing fluxes at the Earth would be respectively 1.2 x 10~^^ and 
2.4 X 10~^^ ergs~^ cm~^, and the detection limits for a point 

ergs~^ cm~^ inafew 



-15 



source for the two instruments are 10" 
lOKsec. 

Recently, Schulz et al. (2003) found that the X-ray emission 
from file MCP-like 07V star 0^ Ori C (Donati et al. 2002) re- 
vealed by Chandra is consistent with a magnetically confined 
wind. Unfortunately, radio observations of this object cannot 
be used to test our model as this star is embedded in the Orion 
nebula, one of the brightest radio sources. 

The capability to reproduce also the X-ray emission repre- 
sents an independent and further confirmation of the validity of 
our model. Multi-wavelength radio observations in conjunction 
with X-ray observations will be a real test for our model. 



5. Conclusion 

The study of radio emission from the MCP stars and, in par- 
ticular, the analysis of the observed modulations offer a unique 
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opportunity to determine the physical properties of the stellar 
magnetosphere. 

The numerical model developed in this paper is used to re- 
produce the 5 GHz radio lightcurves of two well known MCP 
stars: HD 37479 and HD 37017. 

The capability of our numerical model to give us an estima- 
tion of the physical condition of the magnetosphere is limited 
by the lack of multifrequency radio lightcurves. However, the 
spectral information available up to now allows us to focus on a 
small range of possible solutions, all of them indicating a hard 
energetic population of non-thermal-emitting electrons, and an 
inner magnetosphere filled by a thermal plasma consistent with 
a wind-shock model that provides also X-ray emission. 

The possibility to test the radio curves at more than one 
frequency would allow us to put more stringent constraints on 
the model. 

Our simulations provide useful information that may be 
sunomarized as follows: 

we confirm the qualitative model proposed by several 
authors to explain the origin of non-thermal electrons as 
responsible for the observed radio emission from MCP 
stars; 

we point out the importance of the thermal electrons 
trapped in the inner magnetosphere for reproducing the 
rotational modulation of the measured radio flux density 
and the radio spectra of the MCP; this plasma can prroduce 
X-ray emission; 

the length of the current sheets, where the electrons are ac- 
celerated up to relativistic energies, is about one half of the 
Alfven radius; the acceleration process has an efficiency of 
about 10~^, as only this smaU fraction of the electrons in 
the current sheets is accelerated ; those non thermal elec- 
trons have an hard energetic spectrum (A''(7) oc (7 — l)""^), 
with (5 ~ 2; the inner magnetosphere is filled by thermal 
plasma, whose temperature increases outward from 10^ up 
to 10^ K, consistent with a rotating magnetosphere. 

We emphasize the importance of multifrequency radio 
lightcurves and contemporaneous X-ray observations for a fur- 
ther test our model. This will be a powerful investigation tool 
for studying the physics of the stellar magnetosphere. 

Acknowledgements. We thank Carlo Nocita for his help in making 
some figures. We particularly thank the referee Dr. T. Montmerle for 
his constructive criticism which enabled us to strongly improve this 
paper. 

Appendix A: 

A. 1. Procedures 

We numerically developed the physical scenario proposed to 
explain the origin of radio emission arising from the magneto- 
sphere of MCP stars. The aim is to reproduce the centimetric 
radio emission as a function of the rotational phase. 
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For any rotational phase and a given wavelength, the pro- 
cedure is: 

- 3D sampling of the magnetosphere and calculation of the 
magnetic field vector B; 

- definition of the Alfven radius and of inner, middle and 
outer magnetosphere; 

- calculation of the number density rie of the non-thermal 
electrons in each point of the grid; 

- calculation of emission and absorption coefficients; 

- integration of the transfer equation along paths parallel to 
the fine of sight; 

- brightness distribution in the plane of the sky, total flux 
emitted toward the Earth. 

A. 2. The reference frame 

In the oblique rotator model the magnetic axis has an obliquity 
P with respect to the rotational axis 17, which in turn forms an 
angle i with respect to the line of sight. Thus the orientation 
of the magnetosphere changes with respect to the line of sight 
as a function of the rotational phase $. For the purpose of the 
model, the most convenient reference frame where to sample 
the magnetosphere is the fixed frame Ox'y'z', with origin O 
at the center of the star, axis x' direct toward the Earth, axes 
y' and z' in the plane of the sky with z' coinciding with the 
projection of on the plane of the sky. However, since a dipole 
has an axial symmetry, the most convenient reference frame to 
compute the vector B and all the other physical parameters is 
the frame Oxyz, anchored to the star, having axis z coinciding 
with the axis of the dipole. Here the axis x lies in the plane 
defined by the rotational and dipole axes. In this frame, the 
components of vector B are: 

Bx = 3m^ 

By = 3m^ (A.l) 

^- = 4) 

with m the magnetic momentum (m = ^-Bpi?,). If we indi- 
cate by 5ft($) the transformation matrix to pass from the frame 
Oxyz to Ox'y'z', the magnetic field vector in Ox'y'z' is given 
by:' 

B'(r') = 3fiB(5ft- V). (A.2) 

A.3. Rotation of the magnetosphere 

All the physical quantities of the magnetosphere, either vec- 
torial or scalar, are easily computable in the stellar reference 
frame Oxyz. However, we need to know them in the reference 
frame Ox'y'z' of the observer. For this purpose, we need to 
know the transformation matrix of Eq.(A.2). 

The stellar observing system Oxyz is defined by the direc- 
tion of the magnetic axis (z axis), and the intersection between 
the plane of the magnetic equator and the plane where the z 
axis and the rotational axis ft lie (Fig. A. la). 



14 



C. Trigilio et al.: 3D model for radio emission from magnetic CP stars 



In the observing reference frame Ox'y'z' the stellar orien- 
tation is univocally defined by the inclination i of the rotation 
axis with respect to the line of sight, the misalignment /3 of the 
magnetic axis with respect to the rotation axis Q. and the an- 
gle $ corresponding to the rotational phase. The three angles 
i, (3 and $ define the transformation matrix (3, between 
Oxyz and Ox'y'z'. 

The matrix 3? is the result of three space rotations: 



(cos /3 — sin /3 \ 
10 
sin/3 cos/3 / 

rotation of an angle /3 around the y axis (Fig. A.l a-b) so 
that z ^ zi coincides with the rotation axis O; 



(a) (b) 




(e) 



(f) 



cos $ - sin $ ' 
5?2 = sin $ cos $ 
1 

rotation of an angle — <I> around il (Fig. A.l c-d) so that 
xi — > X2 lies in the plane defined by the fine of sight and f2; 



sini cos? 
5R3 = I 10 
— cosi sini 

rotation of {n/2 — i) around y^ axis (Fig. A.l e-f) so that 
X2 — > xz coincides with the Une of sight (axis x'). 

The transformation matrix 3? can so be defined as: 

The inverse transformation, from Ox'y' z' to Oxyz, is then de- 
fined by the three inverse rotations: 

Rotation matrices are unitary, so the inverse rotations are de- 
scribed by the transposed matrices. 

Appendix B: 

B. 1. The gyrosynchrotron emission 

To calculate the total flux density radiated by the magneto- 
sphere, the equation of radiative transfer was numerically inte- 
grated, independently for the ordinary and extraordinary modes 
of propagation, along the axis x'. 

For the gyrosynchrotron emission and absorption coeffi- 
cients rji, and k^,, the approximate expressions given by Klein 
(1987) were used if the harmonic of gyrofrequency was big- 
ger then 4; otherwise the general expressions given by Ramaty 
(1969) have been adopted. 







^={J B 













^ — r 





z' 




B 












Line of sight 







Fig. A.l. Three rotations permit to pass from the stellar reference 
frame Oxyz (panel a) to that of the observer Ox'y'z' (panel f). 

The presence of the star was taken into account by putting 
at the center of the grid a sphere of radius equal to the stellar 
radius (i?*), and opaque to the radiation. 

Defining i,j and k as the index over x', y' and z', the 
method used to inlegrale the equation of radiative transfer may 
be sunmiarized as foUows: 

- calculation of the specific intensity inside each cube ele- 
ment of geometrical depth Al: 



l _ g-fe„(i,j,fe)Ai 



^i'(^) J) k) 

- calculation of the optical depth of the column matter be- 
tween each grid element and the Earth: 



N 



Tu{i,3,k) = ^ k^{i',j,k)Al 



- calculation of the specific intensity 7^ arising from each 
direction parallel to line of sight: 



L{j,k) = ^AI^{i,j,k)e- 



Each element of matrix {j, k) represents the specific intensity 
in the plane of the sky. 

The total flux density can be derived using the relation: 

j k 

where D is the distance of the source. 
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